METHOD FOR REDUCING NOISE, AND 
PULSE PHOTOMETER USING THE METHOD 



BACKGROUND OF THE INVENTION 

The invention relates to a device for measuring physiological 
information through use of photoplethysmograms measured at a plurality of 
wavelengths, and more particularly, to an improvement in a pulse photometer 
used for diagnosing a problem in a respiratory system or a circulatory system. 

A known technique for measuring a concentration of light absorbing 
materials contained in blood during medical care is pulse photometry intended 
for measuring an oxygen saturation Sp02, the concentration of abnormal 
hemoglobin such as carboxyhemoglobin (COHb) and methemoglobin (MetHb), 
or a concentration of injected pigment. A device for measuring oxygen 
saturation Sp02 in particular is called a pulse oximeter. 

The principle of the pulse oximeter is to irradiate onto living tissue 
light rays of a plurality of wavelengths exhibiting different light absorbing 
characteristics with respect to a material of interest and to determine the 
concentration of the material on the basis of a pulse wave produced through 
consecutive measurement of quantity of transmitted light. 

Japanese Patent No. 3270917 (cf., Claims 1, 2 and Figs. 2 and 4) 
discloses determination of a saturation level of oxygen in arterial blood and a 
concentration of light-absorbing materials in the same by: emitting two 
light-rays of different wavelengths onto living tissue; plotting graphs with levels 
of the resultant two pulse waves originating from transmitted light being taken 
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as vertical and horizontal axes, respectively, to thereby determine regression 
lines; and determining the oxygen saturation and the concentration of 
light-absorbing materials in the arterial blood on the basis of gradients of the 
regression lines. This technique enhances measurement accuracy and 
5 reduces power consumption. 

However, a large amount of computing operation is still required to 
determine regression lines and gradients thereof through use of many sampled 
data sets pertaining to pulse waves of respective wavelengths. 

The principle will now be described by taking, as an example, a pulse 
10 oximeter for measuring the oxygen saturation in arterial blood. 

This technique is not limited to a pulse oximeter but is also applicable 
to a device (or a pulse photometer) which, on the basis of the principle of the 
pulse photometry, measures abnormal hemoglobin (carboxyhemoglobin or 
methemoglobin) or a light absorbing material in blood such as pigment injected 
15 in blood. 

Fig. 1 is a block diagram of such a pulse oximeter. The device will 
now be described by taking the pulse oximeter as an example. 

At the time of measurement, a probe 1 is usually attached to an area 
subjected to the measurement. 

20 A processor 8 produces light emission timings for a light emitter 

(LEDs) 2 which produce light rays of different wavelengths, that is, a red-light 
emitting diode (hereinafter abbreviated "R-LED") and an infrared-light-ray 
emitting diode (hereinafter abbreviated "IR-LED"). The R-LED and the 
IR-LED, both being provided in the probe 1, are driven by an LED driver 4, to 

25 thereby alternately emit light. 
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The light emitted by the LEDs 2 passes through the area subjected to 
the measurement to which the probe is attached, and is received by a light 
receiver (photo diode; PD) 3 provided in the probe 1. 

The thus-received signal is converted from light into electricity by the 
5 light receiver 3 and is subsequently converted into a voltage by an input 
section 5. 

Components reflecting optical characteristics of pulsation of the area 
subjected to the measurement appear in this signal in the form of AC 
components. 

10 By a demodulator 6, the signal output from the light receiver (PD) 3 is 

separated and demodulated into a waveform of an infrared light ray (IR) and 
that of a red light ray (R). The waveforms are digitized into signal trains by an 
A/D converter 7. The signal trains are sent to the processor 8, where the 
signal trains are processed, thereby computing, e.g., oxygen saturation SpC>2. 

15 The signal trains that have been digitized by the A/D converter 7 

corresponding to infrared light (IR) and red light (R) form respective measured 
data sets. 

Light absorbance A of the light that has passed through tissue is 
expressed by Equation (1), with the attenuation of light in blood being taken as 
20 Ab and the attenuation of light in other tissue being taken as At. 

A = Ab + At = \n(Iiri) - ln(Iout) = ECD + At ( 1 ) 

wherein: 

25 lout: Intensity of transmitted light 
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lin: Intensity of incident light 
E: Absorption coefficient of blood 
C: Concentration of hemoglobin in blood 
D: Thickness of blood 
5 When the blood thickness has assumed D+AD as a result of pulsation of blood, 
Equation 1 is transformed as follows. 

A + AAb = ln(//w) - ln(Iout) + A ln(/ow/) = EC(D + AD) + At (2) 

10 Subtracting Equation 1 from Equation 2, Equation 3 is obtained. 

AAb = -A \n{Iout) = EC AD (3) 

Equation 3 is defined for wavelengths XI , X2 as follows. 

15 

AAb, = -AlnCAwf,) = E } CAD (4) 
A^6 2 = -Aln(/ow* 2 ) = £ 2 CAD (5) 

Subscripts 1, 2 indicate that the corresponding terms are relevant to 
20 wavelengths X<\ , A*. Provided that Equation 4/Equation 5 = 0>s, we have 

5 A46 2 Aln(/ow/ 2 ) £ 2 ; 

O is determined by measuring Aln(loutl), Aln(lout2), whereby the 
25 oxygen saturation is determined as follows. 
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Sp0 2 = f(O s ) 



(7) 



When the influence of body movement or a like factor other than the 
component reflecting pulsation is superimposed on Aln(loutl), Aln(lout2), light 
5 absorbance of each of the wavelengths is defined as follows. 

AAb[ = -A \n(Iout\ ) = E.CAD + En.CnbDn (8) 
AAb 2 = -A \n{Ioui 2 ) = E 2 CAD + En 2 CnM)n (9) 

10 wherein: 

En: Absorption coefficient of the noise source 

Cn: Concentration of the noise source 

Dn: Change in thickness of the noise source 
A waveform (pulse wave) reflecting a pulsation component is not solely 
15 obtained, and a signal on which noise is superimposed is observed. In this 
case, a light absorbance ratio <E>' is defined as follows: 

0 ._ _ AlnQfoi//;) _ E } CAD + En^CnADn 
" AAb' 2 ~ Mn(Ioui 2 )~ E 2 CAD + En 2 CnADn 

20 Hence, the light absorbance ratio does not coincide with the oxygen saturation 
in arterial blood. 

Fig. 2A shows waveforms which have been subjected to: processing 
in which transmitted light intensity data at respective wavelengths at 
predetermined time intervals are logarithmic-computed; and then a mean of 
25 the logarithmic-computed values is made to zero or high-pass filtering is 
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performed with respect to the logarithmic-computed values. 

Fig. 2B shows a graph through use of waveform data measured 
substantially at the same time when the waveforms shown in Fig. 2A are 
measured (the time is sufficiently shorter than the duration of a frequency 
5 component of a pulse wave) with the amplitude of infrared light being taken as 
a horizontal axis and the amplitude of red light being taken as a vertical axis. 

If the observed data reflect a mere pulse wave component, the graph 
essentially assumes the shape of a straight line. The gradient of the straight 
line shows the light absorbance ratio Os. 
10 However, when noise is superimposed on the data as shown in Fig. 

3A, there is obtained a result into which the light absorbance ratio of noise and 
the light absorbance ratio of pulsation are merged as shown in Fig. 3B, in 
contrast with the case of a pulse wave component. 

As mentioned above, when measured pulse wave data include noise, 
15 an accurate light absorption ratio cannot be measured, and hence noise must 
be reduced. 

Hitherto-known techniques for reducing this noise are a frequency 
analysis technique and an independent component analysis technique. 

However, when the frequency of noise overlaps the fundamental 
20 wave of a signal component of measured pulse wave data or a harmonic wave 
of the same, the frequency analysis technique encounters difficulty in 
extracting a signal component. 

According to the independent component analysis technique, an 
observed signal cannot be dissolved into independent components beyond 
25 measured waveform data. When a plurality of noise sources are present, 
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difficulty is encountered in dissolving the data. 

SUMMARY OF THE INVENTION 

5 It is therefore an object of the invention to provide a method for 

determining a light absorbance ratio by readily reducing noise from an 
observation signal on which noise is superimposed, to thereby extract pulse 
wave data. 

It is also an object of the invention to provide a pulse photometer 
10 using pulse wave data from which noise has been reduced by the above 
method. 

In order to achieve the above objects, according to the invention, 
there is provided a method of processing an observed pulse wave data, 
comprising steps of: 

15 irradiating a living body with a first light beam having a first 

wavelength and a second light beam having a second wavelength which is 

different from the first wavelength; 

converting the first light beam and the second light beam, which have 

been reflected or transmitted from the living body, into a first electric signal 
20 corresponding to the first wavelength and a second electric signal 

corresponding to the second wavelength, as the observed pulse data; 

computing a light absorbance ratio obtained from the first electric 

signal and the second electric signal, for each one of frequency ranges dividing 

an observed frequency band; and 
25 determining that noise is not mixed into the observed pulse wave data 
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in a case where a substantial match exists among light absorbance ratios 
computed for the respective frequency ranges. 

Preferably, the existence of the substantial match of the light 
absorbance ratios is determined with regard to frequency ranges in which at 
least one of the first electric signal and the second electric signal has relatively 
large powers. 

According to the invention, there is also provided a method of 
processing an observed pulse wave data, comprising steps of: 

irradiating a living body with a first light beam having a first 
wavelength and a second light beam having a second wavelength which is 
different from the first wavelength; 

converting the first light beam and the second light beam, which have 
been reflected or transmitted from the living body, into a first electric signal 
corresponding to the first wavelength and a second electric signal 
corresponding to the second wavelength, as the observed pulse data; and 

whitening the first electric signal and the second electric signal by an 
affine transformation using a known light absorbance ratio, in order to separate 
a pulse signal component and a noise component which are contained in the 
observed pulse data. 

Preferably, the affine transformation is performed with the following 

equation: 



where, S is the pulse signal component, N is the noise component, s1 is the 
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first electric signal, s2 is the second electric signal, <|> = Tan" 1 <D, <D is the known 
light absorbance ratio, and 9 is a value selected from a range of 4 to (nl2 - <t>), 
and 

wherein 0 is so selected as to make a norm of the pulse signal 
5 component minimum. 

Here, it is preferable that the signal processing method further 
comprises steps of: 

computing a light absorbance ratio obtained from the first electric 
signal and the second electric signal, for each one of frequency ranges dividing 
1 0 an observed frequency band; and 

determining that noise is not mixed into the observed pulse wave data 
in a case where a substantial match exists among light absorbance ratios 
computed for the respective frequency ranges. 

Here, one of the light absorbance ratios, which are determined that 
15 the noise is not mixed therein, is used as the known light absorbance ratio. 

Preferably, the signal processing method further comprises a step of 
obtaining a signal-to-noise ratio of the observed pulse wave data by performing 
a frequency analysis with respect to the pulse signal component and the noise 
component at each of predetermined frequencies. 
20 Preferably, the signal processing method further comprises a step of 

displaying a pulse wave of the living body, based on the pulse signal 
component. 

Preferably, the signal processing method further comprises a step of 
calculating a pulse rate of the living body based on the pulse signal 
25 component. 

9 



According to the invention, there is also provided a method of 
processing an observed pulse wave data, comprising steps of: 

irradiating a living body with a first light beam having a first 
wavelength and a second light beam having a second wavelength which is 
5 different from the first wavelength; 

converting the first light beam and the second light beam, which have 
been reflected or transmitted from the living body, into a first electric signal 
corresponding to the first wavelength and a second electric signal 
corresponding to the second wavelength, as the observed pulse data; 
10 whitening the first electric signal and the second electric signal to 

separate a pulse signal component and a noise component which are 
contained in the observed pulse data, for each one of frequency ranges 
dividing an observed frequency band. 

Preferably, the step of whitening the first electric signal and the 
15 second electric signal is performed with independent component analysis. 

Preferably, the signal processing method further comprises a step of 
obtaining a signal-to-noise ratio of the observed pulse wave data by performing 
a frequency analysis with respect to the signal component and the noise 
component at each one of the frequency ranges. 
20 According to the invention, there is also provided a pulse photometer 

in which each of the signal processing methods is executed. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The above objects and advantages of the present invention will 
become more apparent by describing in detail preferred exemplary 
embodiments thereof with reference to the accompanying drawings, wherein: 

Fig. 1 is a block diagram showing the configuration of a pulse 
oximeter according to one embodiment of the invention; 

Fig. 2A is a waveform obtained through logarithmeric processing ; 

Fig. 2B is a graph plotted with the amplitude of infrared light being 
taken as a horizontal axis and the amplitude of red light being taken as a 
vertical axis; 

Fig. 3A shows observed (measured) waveforms on which noise is 
superimposed; 

Fig. 3B is a graph plotted with the amplitude of the infrared light being 
taken as a horizontal axis and the amplitude of the red light being taken as the 
vertical axis; 

Fig. 3C is a graph obtained by subjecting the graph shown in Fig. 3B 
to processing performed by the transformation matrix A; 

Fig. 3D is a waveform whose noise has been separated or reduced 
through the processing; 

Fig. 4 is a view showing the relationship between a signal component 
and a noise component; 

Fig. 5 is a view showing an infrared light spectrum and a red light 
spectrum, both obtained through frequency analysis performed with respect to 
the waveforms shown in Fig. 3A; 
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Fig. 6 is a view showing a spectrum obtained through frequency 
analysis of the signal component and the noise component shown in Fig. 3D; 

Fig. 7 shows graphs plotted in the same way as that in Fig. 2B for 
every divided frequency ranges, in a case where noise is superimposed on an 
5 observation signal; 

Fig. 8 shows graphs prepared by the principle shown in Fig. 2B for 
every divided frequency ranges, in a case where noise is not superimposed on 
an observation signal; 

Fig. 9 shows a result obtained by determining a gradient for each of 
10 the frequency ranges shown in Fig. 7, through independent component 
analysis; 

Fig. 10 shows a result obtained by determining a gradient for each of 
the frequency ranges shown in Fig. 8, through the independent component 
analysis; and 

15 Fig. 11 is a flowchart showing a processing flow performed in the 

pulse oximeter shown in Fig. 1. 

DETAILED DESCRIPTION OF THE INVENTION 

20 (How to Separate Signal Component from Noise Component by Utilization of 

Known Light Absorbance Ratio) 

First, the light received by the light receiver 3 provided in the probe 1 

having such a configuration as that shown in Fig. 1 is converted into a voltage 

by the input section 5. 
25 A component reflecting an optical characteristic of pulsation arising in 
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an area subjected to the measurement appears in the received light signal as 
an AC component. 

The signal output from the light receiver 3 is separated and 
demodulated into a waveform of infrared light (IR) and a waveform of red light 
5 (R) by the demodulator 6, and the waveforms are digitized by the A/D 
converter 7. 

A data sequence of the infrared light observation signal IR, a data 
sequence of the red light observation signal R, pulse wave components of the 
observation signal S, and noise components of the observation signal N are 
1 0 respectively taken as: 
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IR = {IRi; i=0, 1,2, ...} 
R = {Ri; 1=0,1,2, ...} 
S = {Si; 1=0,1,2,...} 
N = {Ni; 1=0,1,2, ...} 



(11) 
(12) 
(13) 
(14) 
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As shown in Fig. 4, provided that the gradient of the pulse wave component is 
taken as <}> and the gradient of the noise component is taken as cp, the resultant 
observation signals (IRi, Ri) are expressed as the sum of vectors, i.e., a signal 
component Si and a noise component Ni, as follows: 



7/^ (cos* cos^YV 
^R,) [sin* sin <p\.N 0 



(15) 



Here, defining a matrix W as: 
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w = 



^cos<f> cos^ 



(16) 



^sin^ sin <p j 

If a transformation matrix A which is an inverse matrix of W can be 
caused to act on observation data, the signal component S can be separated 
from the noise component N as represented below. 
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= A 











(17) 



This means that the signal component S is converted into a horizontal axis and 
the noise component N is converted into a vertical axis. The transformation 
matrix A to be used for implementing whitening is determined as: 



L S1 



COS0 COS0> 

sin 
1 



sin <p J 

r 



sin(#>-0) 



- V N 

1/ 

/sin0 J 



sin <p - cos <p 
^-sin^ cos^ j 

f cos^ sin <j>^ 
^-sin^ cos^; 



(18) 



Therefore, we have: 
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r s\ 



1/ 



'tantf 
'sin 0 ) 



f cos0 smflYllP 
v -sin (f> cos^ j^R j 



(19) 



The light absorbance ratio O designates the gradient of a signal 
component. The relationship existing between the gradient <j> and the light 
absorbance ratio O is expressed as <D = tan<|>. 

If a period of time is sufficiently shorter than the duration of a pulse 
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rate, an immediately preceding value (e.g., a value obtained 10 msec, 
previously) can be taken as a known value of <D. Alternatively, the light 
absorbance ratio which has been determined that the noise is not mixed 
therein is used as the known value of O. Still alternatively, a candidate for $ 
5 obtained by the technique described later can also be employed. 
Next, 9 is expressed as: 

0 = <p - 4> (20) 

10 Since cp is unknown, 9is determined as an angle at which a norm in the 
direction of the horizontal axis becomes minimum between to (k/2 - <J>). 

When the observation signal is processed through use of the 
transformation matrix A using the thus-obtained <j>, 0, a signal component 
appears along the horizontal axis, and the noise component appears along the 

15 vertical axis. 

More specifically, the determinant expressed as Equation 19 means a 
coordinate transformation from a first coordinate system in which the observed 
signals IR and R are selected as the coordinate axes shown in Fig. 3B, to a 
second coordinate system in which the signal component and the noise 

20 component are selected as the coordinate axes shown in Fig. 3C. 

Fig. 4 shows observed values (IRi, Ri) at a certain time point t = ti. 
The vector Xi representing the observed values is a synthesis vector of a 
signal component vector Si and a noise component vector Ni. Upon 
performing the above coordinate transformation, the signal component vector 

25 Si is first operated to be made coincident with the horizontal axis of the second 
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coordinate system. Namely, the vector Si shown in Fig. 4 is rotated clockwise 
by an angle of (j>. The value of $ is obtained through the above relationship 
using the known value of <D. 

The coefficient 0 in Equation 19 corresponds to an operation to make 
5 the noise component vector Ni coincident with the vertical axis of the second 
coordinate system. Since the waveform shown in Fig. 3B is a set of the 
observed values, if the values of $ are made constant for the observed values 
at any of the time points, the values of 0 accordingly vary corresponding to the 
respective time points. In order to uniquely determine the value of 0 as the 
10 coefficient, dispersion of values (norm) in the horizontal axis (signal 
component) direction is checked while varying the value of 0 within a range 
from -<J) to (n/2 - The value of 0 is finally determined as a value making the 
norm minimum. 

Fig. 3C shows a graph obtained by subjecting the graph shown in Fig. 
15 3B to processing performed by the transformation matrix A. Further, Fig. 3D 
shows a waveform from which noise has been separated through the 
processing. 

Fig. 5 shows an infrared-ray spectrum (IR) and a red light spectrum 
(R), both being obtained through frequency analysis performed before 
20 processing, and Fig. 6 shows a spectrum obtained through frequency analysis 
of a signal component and a noise component performed after corresponding 
processing. 

(How to Calculate Pulse Rate) 

A peak-to-peak interval Ts (Fig. 3D) is obtained from the signal 
25 waveform extracted through the foregoing conversion processing, whereby a 
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pulse rate = 1/Ts x 60 [beats/min.] can be readily determined. 

Alternatively, the fundamental frequency fs of a pulse wave in relation 
to the spectrum (Fig. 6) obtained from the processed signal through frequency 
analysis becomes clearer than the spectrum (Fig. 5) obtained through 
5 frequency analysis performed before processing. Hence, a pulse rate can be 
computed by a pulse rate = fs x 60 [beats/min.]. 
(How to Determine Whether Noise is Mixed into Observation Signal) 

An observation signal is separated into signals each associated with 
one of predetermined frequency ranges which are divided. The division of the 
1 0 observed frequency band may be performed for consecutive frequency ranges 
or discrete ranges. 

In each frequency range, an infrared light component is taken as a 
horizontal axis (or a vertical axis), and a red light component is taken as the 
axis orthogonal to the horizontal axis (or the vertical axis). 
15 Fig. 7 shows graphs plotted with the infrared light component being 

taken as the horizontal axis and the red light component being taken as the 
vertical axis. 

A gradient corresponding to the principal component of the graph is 
determined over the entire pulse wave frequency range ("all bands" shown in 
20 Figs. 7 and 8) and for each frequency range obtained after separation. The 
observation signal is determined to be a single signal or the sum of a plurality 
of signals, on the basis of whether a match exists between the gradients. 

Specifically, a method for determining whether an observation signal 
is a single signal is to determine whether the gradient of the entire pulse wave 
25 frequency band and the gradients of all the separated frequency range fall 
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within a predetermined range, or whether the gradient of a frequency range 
having a relatively large power and the gradient of the entire pulse wave 
frequency band fall within the predetermined range. 

Specifically, Figs. 7 and 8 show that the entire pulse wave frequency 
5 band and the range from 0.5 Hz to 6 Hz are divided in increments of 0.5 Hz. 
The gradient of the principal component is described at the top of each graph. 

Fig. 7 shows different values in ranges having a relatively large power, 
such as 0.5 to 1 Hz: 0.54, 3.5 to 4.0 Hz: 0.84, and 4.0 to 4.5 Hz: 0.84. 
Summation of a plurality of signals can be determined to be observed. 
10 Fig. 8 shows a value of 0.50 is exhibited in a range of 1 to 1.5 Hz, a 

range of 2 to 2.5 Hz, and a range of 2.5 to 3 Hz, each having a relatively large 
power. Hence, the observation signal is determined to be a single signal. 
(Computation of Candidate for Oxygen Saturation) 

The pulse wave signal frequency band is divided into a plurality of 
15 ranges, as in the case of determination of mixing of noise. Independent 
component analysis is performed for each of the thus-separated frequency 
ranges. 

Under the assumption that a longitudinal vector of the observation 
signal is taken as X and a transposed matrix is taken as T, the following 
20 equation is derived from Equations 1 1 and 12. 



X — ( X 0 , X x , X 2 , • • -) 7 



(21) 



(*,)= 



0 = 0,1,2,-) 



(22) 



25 
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A variance-covariance matrix H of the longitudinal vector X of the observation 
signal is expressed as: 



H = X T -X (23) 

5 

A diagonal matrix having eigenvalues as diagonal elements is taken as A, and 
a matrix formed from eigenvectors is taken as r. Provided that the 
transformation matrix is taken as U, the independent component Y is obtained 

as: 

10 

Y =u-x = r-A 2 -r T -X (24) 

The independent component Y is dissolved into independent components. 

If signal components are present in the thus-separated frequency 
15 ranges, the signal components can be converted into a signal 
component-noise component plane such as that shown in Fig. 3C. 

If the matrix formed from eigenvectors is converted through use of an 
inverse matrix U" 1 of the transformation matrix U, a vector is projected onto the 
IR-R plane. The gradient of the vector expressing the signal component is 
20 obtained as <Ds, and the gradient of the vector representing the noise 
component is obtained as On. 

Fig. 10 shows results of determination of the gradient of the 
observation signal shown in Fig. 2 on each of the frequency ranges shown in 
Fig. 8. The light absorbance ratios 1 are understood to be close to each other 
25 over the respective frequencies (i.e., 0.5+0.02). Fig. 9 shows results obtained 
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by processing the graphs shown in Fig. 7, and a result showing <£>s is obtained 
within a range including fs shown in Fig. 6. The oxygen saturation Sp02 is 
obtained as a function of <I>s; that is, Sp02 = f (<I>s). 

Operation of the processor of the device block shown in Fig. 1 will 
5 now be described in detail with reference to Fig. 1 1 . 

When measurement is commenced, digital signals obtained as a 
result of observation of infrared light and red light are subjected to preliminary 
processing (step S1). Here, the preliminary processing refers to logarithmic 
computation or pulse wave band filtering on the basis of the Lambert-Beer law. 
10 The signals are divided into predetermined frequency ranges (e.g., in 

increments of 0.5 Hz) through use of band-pass filters, thereby determining a 
principal component of light absorbance ratio (i.e., the gradient of the principal 
component on an IR-R plane) for each frequency range (step S2). 

Whether a match exists between the gradients determined for the 
15 respective frequency ranges is determined, thereby determining a candidate 
for a light absorbance ratio of a pulse wave (step S3). A plurality of 
candidates may be acceptable. 

A determination is made as to whether the candidate (signal source) 
is single or plural (step S4), on the basis of the results obtained in the step S2. 
20 In a case where the determination is single, oxygen saturation is 

determined from the principal component light absorbance ratio, and the 
fundamental cycle is determined from the result obtained in the step S1. 
Thus, a pulse rate is computed (step S5). 

In a case where the determination is plural, at the step S6, affine 
25 transformation is performed through use of the candidates for the light 
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absorbance ratio obtained in the step S3 or the result of the light absorbance 
ratios corresponding to an immediately-preceding SpC>2 value (sufficiently 
shorter than the pulse rate: for example, 10 msec faster), or a light absorbance 
ratio which has been determined that the noise is not mixed therein. Equation 
5 19 is used for the transformation matrix, and the matrix is rotated with 
reference to <J> (i.e., a value corresponding to a light absorbance ratio), and 0 is 
determined as an angle at which a norm value in the direction of the <j> axis 
becomes minimum. 

In order to identify a pulse rate, the fundamental cycle Ts (see Fig. 

10 3D) is determined from a cycle determined along the time axis of the signal 
obtained in the direction of $ determined in the step S6, or the frequency 
analysis is performed to identify the fundamental frequency fs (see Fig. 6). A 
signal-to-noise (S/N) ratio is determined for each frequency, and the 
thus-determined ratio is taken as one index for evaluating the reliability of the 

15 obtained light absorbance ratio (step S7). 

On the other hand, time domain data for the respective frequency 
ranges, which are divided at the step S2, are used to perform independent 
component analysis (step S8). 

The gradients at the IR-R plane of the eigenvectors obtained in the 

20 step S8 are determined, and it is checked whether a match exists between the 
frequency ranges, to thus determine a candidate for a light absorbance ratio 
(step S9). A plurality of candidates may be acceptable. 

Outputs obtained in the steps S3, S5, S7, and S9 are compared with 
the light absorbance ratio and bands obtained in the past. Current oxygen 

25 saturation is determined from a common item, consistency, and reliability 



21 



(output from the step S7). The oxygen saturation is averaged, as required, 
and the thus-averaged oxygen saturation is output (step S10). 

When measurement is performed continuously, processing returns to 
the step S1, and processing is iterated. If measurement is not performed 
5 continuously, measurement is terminated. 

According to the invention, the pulse wave signal is separated from 
the noise signal through use of a known light absorbance ratio, thereby 
acquiring a pulse wave signal whose noise has been reduced. 

affine transformation is used for utilization of a known light 
10 absorbance ratio, thereby whitening an observation signal. The pulse wave 
signal is separated from the noise signal. Thus, a pulse wave signal whose 
noise has been reduced can be acquired. 

The pulse wave signal frequency band is divided into a plurality of 
frequency ranges. The principal component of the light absorbance ratio is 
15 determined for each frequency band. By judging whether a match exists 
between the light absorbance ratios, a determination can be made as to 
whether noise is mixed into the observation signal. 

The affine transformation utilizing a known light absorbance ratio 
enables determination of a first rotation utilizing a known light absorbance ratio 
20 and a coefficient suitable for reducing noise. Thus, a pulse wave signal 
whose noise has been reduced can be acquired. 

A pulse can be detected from the pulse wave signal whose noise has 
been reduced, or a pulse rate can be determined through frequency analysis of 
the pulse wave signal. 
25 Dividing the pulse wave signal frequency band into a plurality of 
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frequency ranges enables a reduction in the number of signal sources. 
Further, performing whitening operation on a per-frequency-range basis 
enables determination of a light absorbance ratio of the pulse wave signal and 
a light absorbance ratio of the noise signal. 

Since the pulse wave signal is separated from the noise signal, an 
S/N ratio is obtained for each divided frequency range, and the thus-obtained 
S/N ratio can be utilized as an index to be used for evaluating a signal. 
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